Far-from-equilibrium growth of thin films in a temperature gradient 
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The irreversible growth of thin films under far-from-equilibrium conditions is studied in (2 + 
1)— dimensional strip geometries. Across one of the transverse directions, a temperature gradient 
is applied by thermal baths at fixed temperatures between Ti and T2, where Ti < T^° m < T2 and 
rphom _ q gg(i) j s -j-jjg cr it;ical temperature of the system in contact with an homogeneous thermal 
bath. By using standard finite-size scaling methods, we characterized a continuous order-disorder 
phase transition driven by the thermal bath gradient with critical temperature T c = 0.84(2) and 
critical exponents v = f.53(6), 7 = 2.54(lf), and j3 — 0.26(8), which belong to a different universality 
class from that of films grown in an homogeneous bath. Furthermore, the effects of the temperature 
gradient are analyzed by means of a bond model that captures the growth dynamics. The interplay 
of geometry and thermal bath asymmetries leads to growth bond flux asymmetries and the onset of 
transverse ordering effects that explain qualitatively the shift in the critical temperature. 

PACS numbers: 81.15.Aa, 68.35.Rh, 64.60.De, 05.50.+q 



From an experimental and applied science perspec- 
tive, thin films are extremely useful in a variety of areas, 
from the manufacture of optics (for reflective and anti- 
reflective coatings, self-cleaning glasses, etc.) to electron- 
ics (layers of insulators, semiconductors, and conductors 
from integrated circuits) and packaging (e.g. aluminium- 
coated PET films). Since the growth temperature is one 
of the key parameters in the formation of ordered thin 
films, several experiments have focused on the influence 
of a temperature gradient during film growth. In an 
early experiment by Tanaka et al. [H, magnetic Tb-Fe 
films were grown between two substrates with a tem- 
perature gradient, reporting the observation of perpen- 
dicular magnetic anisotropies and other gradient-driven 
structural features. More recently, Schwickert et al. [2| 
introduced the temperature wedge method where a cal- 
ibrated temperature gradient of several hundred Kelvin 
was established across the substrate during co-deposition 
of Fe and Pt on MgO(OOl) and MgO(llO) substrates. 
Also, Yongxiong et al. [3[ have investigated the evolu- 
tion of Fe oxide nanostructures on GaAs(lOO) by using 
a multi-technique experimental setup. In these studies, 
nanoscale epitaxial Fe films were grown, oxidized, and an- 
nealed using a gradient temperature method, which led 
to nanostripes with uniaxial magnetic anisotropy. More- 
over, the experimental progress on this field has led to a 
variety of technological applications as well [H ■ 

From a theoretical standpoint, the so-called gradient 
percolation method was originally introduced to inves- 
tigate the percolation transition and later applied to a 
variety of problems as e.g. first- and second-order ir- 
reversible phase transitions in far-from-equilibrium sys- 
tems [f|. In magnetic systems, damage spreading pro- 
cesses in a temperature gradient [?J and studies of several 
one-dimensional models Q have been followed by the in- 
vestigation of the kinetic Ising model in two dimensions 



under different kinds of dynamics Q . 

In the context of these recent experimental and theo- 
retical investigations, this work focuses on the irreversible 
growth of magnetic thin films in a transverse tempera- 
ture gradient by means of extensive Monte Carlo simu- 
lations. To our best knowledge, this is the first inves- 
tigation on the growth of magnetic films that presents 
a full characterization of a gradient-driven phase transi- 
tion and universality class. Magnetic films growing un- 
der far-from-equilibrium conditions are investigated by 
using the magnetic Eden model an extension of the 
classical Eden model in which particles have an addi- 
tional degree of freedom (the spin). Notice that, despite 
the fact that we use a magnetic terminology through- 
out, this work is not restricted to magnetic films and can 
be applied to other physical systems, e.g. binary alloys. 
Earlier studies have shown that films growing in (d + 1)- 
dimensional strip geometries in an homogeneous thermal 



bath are noncritical for d = 1 However, for d = 2 

they undergo order-disorder phase transitions that take 
place at T^ om — 0.69(1). Intriguingly, the critical expo- 
nents for this far-from-equilibrium growth model agree 
within error bars with the exact exponents for the equi- 
librium Ising model in d = 2 [ll[ . 

In this work, magnetic films in (2 + 1)— dimensions are 
studied in a square-lattice geometry L x x L y x L z , where 



L z > L x 



L is the growth direction. The starting 



seed for the growing film is a plane of L x L up spins 
(i.e. S — 1) placed at z = 1 and film growth takes 
place along the positive longitudinal direction (z > 2). 
The growth process is characterized by an initial tran- 
sient length It ~ L 2 followed by a nonequilibrium sta- 
tionary state that is independent of the initial configu- 
ration [ll|. Therefore, any choice for the seed leads to 
the same stationary states for z ^> By disregarding 
the transient region, all results reported in this paper are 
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obtained under stationary conditions. Across one of the 
transverse directions (the y— axis), a temperature gradi- 
ent is applied by thermal baths at fixed temperatures 
linearly varying between Ti and T 2 . Hence, we adopt 
open boundary conditions along the y— direction. Across 
the other transverse direction (the a;— axis), continuous 
boundary conditions are considered, meaning that sites 
at (x = l,y,z) and (x = L,y,z) are lattice neighbors. 
Notice that the assumption of continuous boundary con- 
ditions is the standard method used in Monte Carlo sim- 
ulations to eliminate the impact of surface effects that 
would arise otherwise due to missing neighbors in the 
lattice boundaries [l2j. Unless stated otherwise, results 
below correspond to T\ = 0.5 and T 2 = 1.5. We checked 
that changing the gradient has only mild finite-size ef- 
fects that become negligible in the thermodynamic limit. 
That is, by applying different gradients the system's ob- 
servables such as the magnetization, susceptibility and 
higher-order cumulants are shifted in the transition re- 
gion by small amounts, which vanish in the L — > oo limit. 
More details with a full account of finite-size effects will 
be published elsewhere 13 1. 

Films are grown by selectively adding spins (S r = ±1) 
to perimeter sites, which are defined as the nearest- 
neighbor (NN) empty sites of the already occupied ones. 
Let us recall that the substrate is a 3D cubic lattice and 
therefore each lattice site has 6 NN sites. The deposi- 
tion of new spins is irreversible (i.e. once added, spins 
do not flip, detach, nor diffuse). Considering a ferromag- 
netic interaction of strength J > between NN spins, 
the energy of a given configuration of spins is given by 
E = — (J/2) J2( r T n S r S r >, where the summation is taken 
over occupied NN sites. The Boltzmann constant is set 
equal to unity throughout, and both temperature and 
energy are measured in units of J. The probability for 
a perimeter site to be occupied is proportional to the 
Boltzmann factor exp(— AE/T), where AE is the change 
of energy involved in the addition of the new particle and 
T is the temperature at the perimeter site. At each step, 
all perimeter sites have to be considered and the proba- 
bilities of adding a new (either up or down) spin to each 
site must be evaluated. Since this procedure requires up- 
dating the probabilities at each time step, the algorithm 
is very slow compared with those used in equilibrium 
spin models. Involving a significant computational ef- 
fort, clusters having up to 10 9 spins have been generated 
for lattice sizes in the range 12 < L < 96. As in the case 
of the classical Eden model, the magnetic Eden model 
leads to a compact bulk and a sclf-affinc growth inter- 
face (loj . Snapshots and further details on the growth 
behavior will be given elsewhere (l3j |. 

In order to study magnetic films growing in a gradi- 
ent, the appropriate order parameter is the mean ab- 
solute magnetization of transverse columns at constant 
temperature, i.e. (H)(y) = (1/L) (| £) x S xyz \), where 
(...) denotes averages along the growth direction z, and 




FIG. 1: Binder cumulant as a function of the layer tem- 
perature for different system sizes, as indicated. The in- 
set shows the cumulant intersections for the larger systems 
(32 < L < 96), which determine T c = 0.84(2). 



the absolute value avoids shortcomings due to finite-size 
effects, as in standard Monte Carlo simulations [l2T ]. 

The Binder cumulant, Ui(y) = 1 — (m 4 ) /3(m 2 ) 2 , is 
useful to locate the critical temperature because, in the 
low-temperature ordered region, C/4 tends to the value 
2/3, while in the high-temperature disordered region, U4 
tends to 0. Thus, in the thermodynamic limit, U4 be- 
comes discontinuous exactly at the critical temperature 
14j . Figure [1] shows U4 as a function of the layer temper- 
ature for different system sizes. Notice that, by our setup 
assumptions, each layer at fixed y is subjected to a con- 
stant temperature T{y) =T t + (T 2 - T t )(y - 1)/(L - 1) 
maintained by a thermal bath. We assume that each 
layer's temperature can be maintained at the same value 
through the film's growth process (that is, we assume 
that the deposition process does not affect the layer's 
temperature, which is fixed by that layer's thermal bath). 
The results from Figure Q] show evidence for the existence 
of a continuous phase transition driven by the tempera- 
ture gradient. Indeed, if a non-trivial critical tempera- 
ture exists in the thermodynamic limit, we should expect 
the cumulants for different lattice sizes to intersect near 
the critical temperature 14{. The inset in Figure [Tj shows 



a detailed view of the data for the largest lattice sizes 
available, where the intersection region is also indicated. 
Based on this observation, we determine the critical tem- 
perature in the thermodynamic limit as T c = 0.84(2). 

In the following, we apply standard finite-size scal- 
ing techniques 12, [lj| to determine the critical expo- 
nents that characterize the system's critical behavior 
and universality class. Figure [5] shows a log-log plot 
of the finite-size pseudo-critical temperature T C (L) as a 
function of the inverse of the system linear size, where 
T C (L) is defined as the temperature corresponding to 
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FIG. 2: (Color online) Log-log plot of the finite-size pseudo- 
critical temperature T C (L) as a function of the inverse of the 
system linear size, L . Finite-size scaling fits to the data 
obtained by using T c = 0.82, 0.84, and 0.86 lead to the critical 
exponent v = 1.53(6). Inset: plot of (\m\) xL w " vs \T-T C \ x 
L 1 '" (with ft = 0.26) showing a data collapse for different 
system sizes in the range 32 < L < 96. 

(|m|) = 0.5. The finite-size scaling theory predicts that 
\T C — T C (L)\ oc L^ 1 /", where v is the exponent that char- 
acterizes the divergence of the correlation length at crit- 
icality. Least-squares fits to the data obtained by using 
T c = 0.82, 0.84, and 0.86 are also shown. We observe 
that the finite-size scaling relationship fits the data very 
well within the range of values for the critical tempera- 
ture that was derived from the intersection of Binder's 
cumulants. From these fits, we obtain the critical expo- 
nent v — 1.53(6), where the error bars reflect the error 
derived from the evaluation of T c as well as the statistical 
error. On the other hand, the finite-size scaling theory 
predicts that plots of {\m\)L^/ v vs \T-T C \L X I V for differ- 
ent lattice sizes should collapse near the critical region. 
The inset to Figure [5] shows the data collapse obtained by 
using (3 = 0.26 (that is determined from the hyperscaling 
relation, as explained below) with two separate branches 
corresponding to the low- and high-temperature regions. 

Let us now calculate the critical exponent 7, which de- 
scribes the divergence of the susceptibility at the critical 
point. The susceptibility is defined as \ = (i 2 /T) x 
((to 2 ) — (|m|) 2 ). From the finite-size scaling theory (15| . 
the exponent ratio ^/v is related to the peak of the 
susceptibility measured in finite samples of size L by 
Xmax oc L 1 ^ . Figure [3] shows the maxima of \ vs 
L, where the solid line is a fit to the data that yields 
7/1/ = 1.66(3), where the error bars reflect the statisti- 
cal error from the fit. By using this ratio and the value 
already obtained for v, we determine 7 = 2.54(11). The 
insets to Figure display plots of xL~ 1,u vs \T -T C \L X I V 
for different lattice sizes in the range 32 < L < 96. Using 



FIG. 3: (Color online) Log-log plot of the susceptibility max- 
ima as a function of the system linear size, where the solid 
line is a finite-size scaling fit that yields y/u = 1.66(3). The 
insets display plots of \L vs \T — T C \L X I V showing sep- 
arately the data collapse for (a) the low-temperature branch 
and (b) the high-temperature branch. 

the critical temperature as determined by the susceptibil- 
ity peaks, the data collapse is shown separately for (a) the 
low-temperature branch and (b) the high-temperature 
branch. In the former case, data from low-temperature 
layers near T\ = 0.5 depart from the collapse and have 
been removed. However, the collapse near the critical 
region is remarkable and agrees very well with the ex- 
pectations from the finite-size scaling theory. 

By replacing the exponents v and 7 in the hyperscaling 
relation dv — 2(3 — 7 = with d = 2, we determine the 
exponent j3 = 0.26(8). Indeed, we anticipated this value 
of f3 when we considered the data collapse of the scaled 
magnetization (see the inset to Figure [5] above) . The 
excellent data collapse near the critical region confirms 
the consistency and robustness of the obtained results. 
Based on these findings, we conjecture that magnetic 
Eden films grown in a temperature gradient belong to 
a new universality class characterized by critical expo- 
nents v = 3/2,7 = 5/2, and f3 = 1/4. In contrast, previ- 
ous studies found that the critical exponents for magnetic 
Eden films grown in an homogeneous bath agree within 
error bars with the exact exponents for the Ising model 
in d = 2 [IH, namely v = 1,7 = 7/4, and (3 = 1/8. 

Let us now explore the growth dynamics by means 
of a bond representation. For this purpose, to each 
pair of neighboring sites we assign a directed bond that 
points from the earlier occupied site to the later occupied 
one. The components of the bond flux field </> at a site 
s = xi + yj + zk are defined as 4> x (s) = b[s, s + i], 4>y{s) = 
b[s, s+j], and z (s) = b[s,s + k], where b[s{,S2\ = +1 
(— 1) if the bond points from s{ to sa (sa to s{). Figure^] 
shows the components of the mean bond flux as a func- 
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FIG. 4: Mean growth bond flux components as a function 
of the gradient span AT for L = 32 and Ti = 0.5. Asym- 
metries due to the temperature gradient and the substrate 
geometry are responsible for net bond fluxes along the y and 
z directions. Along the tranverse direction x the system is 
fully symmetric, so no net bond fluxes are observed. 

tion of the gradient span AT = Ti — T\ for L = 32 and 
Ti = 0.5. As expected from the symmetry along x, one 
has that ((j) x ) = regardless of AT. For AT = 0, the 
system is also symmetric along y, so no net bond flux is 
observed. When a gradient is applied, however, this sym- 
metry is broken. Since the growth probabilities depend 
on the Boltzmann factor exp(— AE/T(y)), where T(y) is 
the layer's temperature, the thermal asymmetries intro- 
duced by the gradient favor spin deposition on the colder 
layers. This phenomenon is captured by the observed 
net bond flux {(f) y ) > 0. Indeed, as shown in Figure IH 
the thermal asymmetries cause (<j) y ) to grow steeply up 
to ((f>y) m ax ~ 0.75 followed by a moderate decrease for 
larger gradients, which is due to the onset of bulk disorder 
within the hotter layers. Since the net transverse growth 
bond flux is directed from the ordered (cold) layers to- 
wards the disordered (hot) ones, this gradient-induced 
transverse ordering mechanism causes the system's crit- 
ical temperature to increase from T^ om — 0.69(1) to 
T c = 0.84(2). On the other hand, for AT = 0, (<p z ) = 1 
due to the longitudinal asymmetries in the substrate (i.e., 
the semi-infinite strip geometry constrains the system to 
grow along the z > direction). However, when the 
transverse gradient is applied, two effects contribute to 
decrease {<j> z ): (i) the onset of the transverse bond flux, 
which creates transverse domains in the active perimeter 
and causes some of the added spins to grow backwards; 
(ii) the bulk disorder induced in the hotter layers (which 
also causes {<j> y ) to decrease, as discussed above). 

As a summary, we studied the growth of thin films un- 
der far-from-equilibrium conditions in (2+l)-dimensional 
strip geometries, where a temperature gradient is ap- 



plied across one of the transverse directions. In order 
to investigate the system's critical behavior, we applied 
the finite-size scaling theory. The critical temperature is 
T c = 0.84(2) and the critical exponents are v = 1.53(6), 
7 = 2.54(11), and /3 = 0.26(8). Based on these find- 
ings, we conjecture that the exact exponents are v = 3/2, 
7 = 5/2, and (3 — 1/4, which satisfy the hyperscaling re- 
lationship and point out the existence of a new universal- 
ity class for film growth in a thermal gradient. Moreover, 
we investigated the system's growth dynamics by means 
of a bond model. We found that the interplay of geome- 
try and thermal bath asymmetries leads to growth bond 
flux asymmetries and the onset of transverse ordering ef- 
fects that explain qualitatively the shift observed in the 
critical temperature. 

We believe that this work provides the first full charac- 
terization of a gradient-driven phase transition of growing 
films under far-from-equilibrium conditions. Given the 
great experimental and theoretical interest in thin films 
growing in temperature gradients, as well as the wide va- 
riety of technological applications that benefit from these 
efforts, we hope that this paper will contribute to the 
progress of this research field and stimulate further work. 
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